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ABSTRACT 

We have used a numerical simulation of a turbulent cloud to synthesize maps 
of the thermal emission from dust at a variety of far-IR and sub-mm wave- 
lengths. The average column density and external radiation field in the simula- 
tion is well matched to clouds such as Perseus and Ophiuchus. We use pairs of 
single-wavelength emission maps to derive the dust color temperature and col- 
umn density, and we compare the derived column densities with the true column 
density. We demonstrate that longer wavelength emission maps yield less bi- 
ased estimates of column density than maps made towards the peak of the dust 
emission spectrum. We compare the scatter in the derived column density with 
the observed scatter in Perseus and Ophiuchus. We find that while in Perseus 
all of the observed scatter in the emission-derived versus the extinction-derived 
column density can be attributed to the flawed assumption of isothermal dust 
along each line of sight, in Ophiuchus there is additional scatter above what can 
be explained by the isothermal assumption. Our results imply that variations 
in dust emission properties within a molecular cloud are not necessarily a major 
source of uncertainty in column density measurements. 

Subject headings: ISM: clouds — dust, extinction — surveys 



1. Introduction 

In the absence of direct meausurements of extinction, the column density of dust in 
molecular clouds is often traced by dust emission. In many cases the data are maps at two 
wavelengths in the FIR or sub-mm, and the column density is derived by assuming that dust 
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along each line of sight is isothermal (Wood, Myers & Daugherty 1994; Schlegel, Finkbeiner 
& Davis 1998; Schnee et al. 2005). However, each line of sight through a molecular cloud 
is likely to intersect regions with various amounts of shielding, and therefore a range of 
temperatures, so the isothermal assumption may lead to significant errors in the calculation 
of column density. 

In this paper we use a simulation of an externally heated turbulent molecular cloud to 
produce maps of column density and thermal emission at various FIR and sub-mm wave- 
lengths. We discuss the bias in the emission- derived column density as a function of ex- 
tinction and wavelength, and we compare the scatter in the emission-based column density 
estimate with that observed in Schnee et al. (2005) in the Perseus and Ophiuchus molecu- 
lar clouds. We find that assuming isothermality can account for all of the observed point 
by point scatter in comparisons of extinction and emission-based column density maps of 
Perseus without the need to invoke variations in the dust properties within a molecular cloud 
or variations in the emission from stochastically heated grains. In comparison, the scatter in 
Ophiuchus is significantly larger than can be accounted for by the error introduced by the 
isothermal dust assumption alone. 

2. The Model Cloud 

The method of determining the dust distribution, radiative transfer and dust heating, 
only briefly described in this Section, are described in more detail in Bethell et al. (2004). 

2.1. Density Structure 

The density structure is created by a three-dimensional simulation of driven MHD tur- 
bulence on a cubical grid with 128 cells on a side (Heitsch et al. 2001). Although the original 
ensemble of simulations was created with widely varying physical parameters (Mach number, 
plasma beta etc.), the radiative transfer through the resulting structures is largely insensi- 
tive to these differences. To simulate Perseus and Ophiucus we use one such model and 
scale the density and size of our model cloud so that the column density map has an average 
Ay ~ 2, which is the median Ay in our Perseus and Ophiuchus extinction maps. This choice 
of model is meant to be illustrative of the broad effects of including 'realistic' dumpiness 
in radiative transfer calculations of molecular clouds. When compared to a uniform density 
distribution of similar mass the inclusion of dumpiness tends to expose more of the volume 
to the interstellar radiation field (ISRF) while on average the mass is both less illuminated 
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and sees bluer light. In the clumpy case there is also little correlation between the intra- 
cloud radiation field at a point and its proximity to the cloud edge. Instead, there is a large 
point-to-point scatter, arising from the details of the dumpiness (Bethell et al. 2004). 

2.2. Radiative transfer 

A Reverse Monte Carlo radiative transfer code (Bethell et al. 2004) is used to calculate 
the penetration of an external isotropic and homogeneous ISRF, here the ISRF given by 
Mathis, Mezger & Panagia (1983). The reverse Monte Carlo technique enables the user to 
control the Monte Carlo photon noise in regions which are well shielded, in particular the 
clumps where most of the mass resides. The re-emission by dust of sub-mm radiation is 
assumed to be optically thin throughout our cloud. 

2.3. Dust properties 

To calculate self consistently the attenuation of the ISRF by dust, the subsequent equi- 
librium dust temperature and the FIR emission we use astronomical silicate and graphite 
grain opacities (Draine & Lee 1984) 1 and a power-law grain size distribution (Mathis, 
Rumpl & Nordsieck 1977). These are chosen to reproduce an extinction curve with Ry = 
Ay / E(B — V) = 3.1. No attempt is made to include very small grains which are heated tran- 
siently to high temperatures, nor do we include polycyclic aromatic hydrocarbon emission. 
Both these populations contribute significantly to wavelengths shorter than 100 /xm. 

At far infrared wavelengths (A > 100/zm) the dust ensemble emission is described by a 
blackbody spectrum modified by a power law (S\ oc B\(T d ) A _|8 ) with an emissivity spectral 
index of (5 — 2. With the dust ensemble fully prescribed we can calculate the equilibrium 
grain temperatures of each grain type and grain size at each location. Not only do dust 
temperatures vary along each line of sight, they also vary amongst grain types and sizes at 
any point within the cloud. 



For tabulated optical properties visit the website www.astro.princeton.edu/~draine. 
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3. Derived Color Temperature and Column Density 

The method used to determine the dust temperature and column density from two 
emission maps is similar to that used in Schnee et al. (2005). The dust temperature (T^) is 
determined by the ratio of two flux densities and the column density can be derived from 
either flux and the derived color temperature of the dust. The pairs of fluxes that we study 
in this paper are 60 and 100 /im (the longest wavelength IRAS bands), 100 and 240 /im 
(those used in Schlegel, Finkbeiner & Davis (1998), hereafter SFD98), 100 and 850 /im and 
450 and 850 /im (wavelengths observed by SCUBA on the JCMT (Holland et al. 1999)). 

To determine the dust color temperature, we use the equation: 

= ( AA ~ (5+/3) f exp[hc/X 2 kT d ]-l \ 

\\J \exp[hc/ ■A 1 fcT <l ] - 1 J U 
where R is the ratio of the fluxes and (3 is the emissivity spectral index of the dust. We 
calculate the dust temperature with the assumption that (3 = 2, which is true of the dust in 
our model. Note that in using the ratio of two flux maps to get the color temperature, we 
are assuming that a single temperature, T d , can be used to describe all the dust along each 
line of sight. 

The optical depth is determined by the ratio of the flux to the Planck function: 

with the assumption that the dust emission is optically thin. 

The optical depth can then be converted to V-band extinction using: 

Ay = Xr x (3) 

where X is a parameter relating the thermal emission properties of dust to its optical ab- 
sorption qualities. The extinction, Ay, is derived from the column density of dust in our 
model, using the conversions N H /E(B — V) = 5.8 x 10 21 cm~ 2 and A V /E(B — V) = 3.1. We 
use the IDL procedure AMOEBA to determine the value of X that minimizes the difference 
between the emission- derived column density and the true column density in the simulation. 



4. Results 
4.1. Dependence on Wavelength 



The column density of dust derived from pairs of thermal emission maps depends 
strongly on the two wavelengths considered. Figure 1 shows that the 60 and 100 /im derived 
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column density is a linear tracer of the actual column density for Ay < 1 mag, but system- 
atically underestimates the true column density for Ay > 1 mag. The same trend is seen 
in the column density derived from the 100 and 240 /iin emission maps, though the cutoff 
extinction above which the column densities are underestimated is Ay ~4 mag. The column 
density derived from the 100 and 850 /xm maps traces the true column density all the way 
to 10 mag Ay, as does the column density derived from the 450 and 850 /im emission maps. 
The errors introduced by the assumption of isothermal dust along each line of sight for each 
pair of maps is shown in Figure 1. 

The underestimation of column density at Ay > 1 as derived from the 60 and 100 /im 
emission maps comes from an overestimate of dust temperature along lines of sight that pass 
through various environments with different temperatures. The emission is then dominated 
by the warmer regions, while the column density is dominated by the cooler and denser 
regions. The derived color temperature is thus more representative of the warmer regions, 
and the derived optical depth is therefore too low. This effect can be seen clearly in the 
leftmost panel of Figure 1. Because X has been scaled to minimize the difference between 
the emission and absorption derived column densities, the derived column densities of points 
below Ay = 1 are overestimated. This problem is much less significant when the emission 
maps are in the Rayleigh- Jeans portion of dust emission spectrum (Bethell et al. 2004), which 
explains why the two panels on the right side of Figure 1 show a much tighter correlation 
between the derived column density and the true column density than the two panels on 
the left. Similarly, the longer wavelength emission maps, shown in Figure 2, trace the true 
column density much better than the shorter wavelength maps. 



4.2. Variable Dust Properties 

In order to isolate the effect of the assumption of isothermal dust along each line of sight 
through a molecular cloud on the emission- derived column density, we have kept the dust 
properties in our model invariant. However, recent observations have shown that the overall 
submillimeter emissivity of dust grains varies within a molecular cloud (Stepnik et al. 2003) 
and that the emissivity spectral index of dust within a molecular cloud also varies (Dupac et 
al. 2003). The scatter between emission-derived column density and the true column density 
shown in Figure 1 arises solely from the assumption of isothermal dust, and does not include 
the effects of variable dust emission properties. When the column density is derived from 
60 and 100 /im observations, the effect of variable grain properties may not be important, 
because the assumption of isothermal dust makes deriving the column density from emission 
maps impractical above Ay > 1 while dust properties are not seen to vary significantly below 
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Ay < 2 (Stepnik et al. 2003). When the emission maps are longer wavelength (e.g. from 
DIRBE, PRONAOS or SCUBA), the assumption of isothermal dust is less problematic (as 
shown in Figure 1) and the variable dust properties may become the dominant source of 
error in the emission- derived column density. 

4.3. Emission from Very Small Grains 

We have shown in Schnee et al. (2005) that the 60 /im emission of molecular clouds 
comes predominantly from transiently heated very small dust grains (VSG's). In order to 
derive the dust temperature, we have removed the VSG contribution to the 60 /im IRAS 
flux. Because we have no way to determine the VSG contribution on a pixel-by-pixel basis, 
we assume that the percent emission from VSG's at 60 /im is constant everywhere in Perseus 
and Ophiuchus. We then determine the 60 /im scale factor such that the median temperature 
in our maps is equal to the median temperature derived by SFD98 for the same area. The 
SFD98 temperature maps were derived using 100 and 240 /im data and are therefore unbiased 
by the emission from the VSG's. This procedure is identical to that performed in Schnee 
et al. (2005). The 60 /im flux from our model was not subjected to such a scaling factor 
because no VSG contribution to the flux was included. 

Though scaling the 60 /im flux by a constant multiplicative factor can give the proper 
average temperature, each line of sight through our observed molecular clouds should ideally 
be treated with its own scaling factor. The variable VSG contribution to the 60 /im flux 
can't be removed on a pixel-by-pixel basis because the SFD98 resolution is too coarse, which 
introduces errors in the derivation of dust temperature and column density. Applying a 
constant scale factor to the 60 /im flux does not correct for this error, nor does it introduce 
additional scatter in Figure 3. The derived value of X does depend on this scale factor. 

If high resolution data at 100 /tin and longer wavelengths were available for a molecular 
cloud, then one could use those maps to avoid having to subtract the VSG contribution to 
the observed emission. Unfortunately, the dust temperature and column maps produced by 
SFD98 using existing 100 and 240 /im data are of much lower resolution (~1 degree) than 
the IRAS maps (~5'). Even nearby molecular clouds require the higher resolution provided 
by IRAS, so we compare our model molecular cloud with the observations of Perseus and 
Ophiuchus using IRAS emission maps. 
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4.4. Comparison with Observations 

In real molecular clouds the best estimate of column density often comes from the 
measured NIR colors of background stars and is calculated from the reddening done by 
the intervening dust. Column density maps made this way have a fundamental uncertainty 
because the background stars have a spread in their intrinsic NIR colors, and because the 
stars included in each pixel of the column density map sample slightly different lines of sight 
through the molecular cloud. The level of uncertainty is approximately 0.2 mag Ay, and it 
is independent of column density (Ridge et al. 2005). 

In order to compare our simulation to observations, we haved added normally-distributed 
noise (with la = 0.2 mag) to the simulated column density. Assuming (3 = 2, we then de- 
rived the value of X that best matches the emission- derived column density (from the 60 
and 100 /im maps) to the new "true" column density. The result is shown in Figure 3. 

The left panel of Figure 3 appears similar to the Perseus and Ophiuchus extinction 
scatter plots, produced as described in Schnee et al. (2005), shown here in the center and 
right panels. The emission-based column densities in Perseus and Ophiuchus are derived from 
IRAS 60 and 100 jum maps and the extinction is derived from from NIR color excesses (Ridge 
et al. 2005). The column derived from FIR emission maps is overestimated at low extinctions 
and underestimated at high extinctions in both this paper and in observations. The la width 
of the gaussian fit to the ratio of emission-derived column density and absorption-derived 
column density are 0.3 in our simulation, 0.3 in Perseus and 0.4 in Ophiuchus, as shown in 
Figure 3. 

We are able to reproduce the scatter observed in Perseus between the column density 
derived from 60 and 100 /im emission maps and the extinction in our simulation of a turbulent 
molecular cloud. The dust in our simulation does not have variable emission properties, nor 
is there a population of transiently heated VSG's. This suggests that any variations in dust 
properties within Perseus do not dominate the uncertainty in the emission-derived column 
density, nor does the emission from the VSG's. Instead, we find that in Perseus the scatter 
is primarily explained by the assumption of isothermal dust along each line of sight. In 
contrast, the Ophiuchus molecular cloud has a larger emission/absorption derived scatter 
than can be explained solely by the errors introduced by the assumption of isothermal dust 
along each line of sight. 

We have shown in Schnee et al. (2005) that in order to produce column density maps 
of nearby molecular clouds from IRAS emission maps, it is necessary to determine both the 
fraction of 60 /im emission that comes from transiently heated dust grains and the conversion 
factor (A) from 100 /im optical depth to visual extinction, because both of these quantities 
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vary from place to place. These parameters can be determined by the temperature maps of 
SFD98 and 2M ASS/NICER extinction maps. 

5. Summary 

We have shown that column density estimates from longer wavelength thermal emission 
maps are less biased and have less uncertainty than estimates made from shorter wavelength 
maps. In particular, maps of column density derived from 60 and 100 /im emission maps 
or 100 and 240 /zm emission maps systematically overestimate the dust temperature at high 
extinctions (>1 mag Ay and 4 mag A v , respectively) and underestimate the column density. 

We have also shown that the observed uncertainty in column density derived from 
IRAS emission maps of Perseus can be fully accounted for by errors introduced by the 
assumption of isothermal dust along each line of sight through the cloud and uncertainty 
in the extinction-derived column density. We find no evidence of variable dust emission 
properties or variations in the fraction of emission at 60 fj,m from transiently heated VSG's 
creating significant scatter in column density estimates in Perseus. The observed scatter in 
the emission/absorption derived column density in Ophiuchus can be largely, but not fully 
explained by the assumption of isothermal dust along each line of sight, suggesting that the 
VSG emission and/or dust properties vary within the Ophiuchus molecular cloud. 
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Fig. 1. — The column density derived from two emission maps plotted against the true 
column density. The wavelengths used are, from left to right: 60 & 100 /im, 100 & 240 
/im, 100 & 850 fim and 450 & 850 /mi. The values of X that convert from the FIR optical 
depth to visual extinction are r wo = 775, t wo = 231, t wo = 69, and r 45 o = 681. The values 
of n and a refer to the center and width of the gaussian fit to the distribution of the ratio 
(A v . . I Ay, ). 

V v emission I *true / 
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Fig. 2. — The flux maps at the five wavelengths considered in the paper, plotted along with 
the true column density. All of the maps are scaled to have the same dynamic range with 
their maxima black and minima white. 
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Fig. 3. — The column density derived from 60 and 100 fim emission maps (assuming (3 = 2) 
plotted against the "true" column density. The model's true column density has been altered 
by gaussian "noise" with a la of 0.2 mag Ay to match the noise in observational estimates of 
extinction. The middle and right panels come from Schnee et al. (2005), where the extinction 
is measured from NIR colors of backgroud stars and the 60 and 100 /im flux from the IRAS 
survey are used to estimate the column density from dust emission. The blue (dashed) 
line shows the best fit between the emission- derived column density and the "true" column 
density. The slope of the model fit is 0.53 which is intermediate between the Perseus and 
Ophiuchus slopes (0.66 and 0.46, respectively). The values of \i and a refer to the center 
and width of the gaussian fit to the distribution of the ratio (Ay I Ay. .. ,. ). 

O \ v emission I y true .extinction / 



